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At low temperatures, spin dynamics in ideal spin ice is due mainly to dilute, thermally excited magnetic 
monopole excitations. I consider how these will affect the dynamics of a nuclear spin (the same theory applies 
to muon spin resonance if implanted muons do not diffuse). Up to the time scale for nearby monopoles to 
be rearranged, a stretched-exponential form of the relaxation functions is expected. I work out the expected 
exponent in that exponential and the formulas for the Ti (longitudinal) and T2 (dephasing) relaxations, as a 
function of the monopole density (and implicitly the temperature). 
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Introduction — 

The dipolar spin ice state [1-5] occurs at sufficiently low 
temperatures in certain highly frustrated pyrochlore structure 
magnets such as Dy2Ti207. The magnetic ions (here Dy+'^) 
sit on a a "pyrochlore" magnetic lattice consisting of corner- 
sharing tetrahedra, and have a strong uniaxial anisotropy 
along their local (111) symmetry axis. The practically degen- 
erate lowest energy states consist of the (exponentially many) 
configurations in which every tetrahedron has two spins point- 
ing inwards and two spins pointing outwards. The elemen- 
tary excitations in such a material are defect tetrahedra with 
three spins in and one out, or the reverse, which were shown 
to behave as (emergent) magnetic monopole f3\. The low- 
temperature thermodynamics and dynamics of dipolar spin ice 
depend on the density and mobility of the monopoles. 

Measurements of the specific heat and (static) spin cor- 
relations in the presence of a large field |7J, corroborate the 
monopole picture. An analog of the Wien effect in elec- 
trolytes was invoked to infer the effective monopole charge 
from muon spin relaxation |@] but this has been called into 
question |l 9 l ll0tl . The dynamics in magnetic relaxation experi- 
ments pl'l shows an activated dynamics, but not the activation 
energy expected from the basic monopole model [12]. Low- 
temperature nuclear magnetic resonance (NMR) jisl il4tl or 
muon spin resonance experiments | 15 , T^l) may aid in dis- 
entangling this confusion, as they probe local spin dynamics 
rather than the uniform bulk magnetization. 

In this Letter, I will first review the essentials of spin ice and 
its emergent monopoles, and then set up the simplest possible 
model for NMR relaxation functions, both the Ti (longitudi- 
nal) and T2 (dephasing) kind. In the presence of disorder that 
is fixed for the duration of the measurement — which includes 
monopoles if they do not diffuse too fast — the NMR signal 
is a superposition of relaxation functions due to the inequiv- 
alent environments of many probe spins, giving a stretched 
exponential form. In the monopole case, the form of the Ti 
relaxation function will be shown to have the same stretched- 
exponential form as the case of fixed magnetic impurities; but 
the T2 relaxation function is novel, in that a power of t^/^ is 
replaced by t'^/^ in the exponential. Experiment does find 
stretched exponentials, but I will show their parameters are 
incompatible with a monopole origin. It is proposed that in- 
terstitial "stuffed" moments (additional to the spin ice Hamil- 



tonian) are responsible for the NMR signal at very low tem- 
peratures. 

Dipolar spin ice and dumbbell model 

The experiments in Ref. [Ts^ are on Dy2Ti207 in which the 
magnetic ions are rare earth Dy^'^ (Ho2Ti207 is similar). 

Each Dy spin has a strong Ising-Iike anisotropy with easy 
axis aligned along the (111) direction of its local 3 -fold sym- 
metry. The pyrochlore lattice nearest-neighbor distance is 
flp = 1/ a/S times the cubic lattice constant. The tetrahedron 
centers, separated by ad — ■\/3/4 times the cubic lattice con- 
stant, constitute the vertices of a diamond lattice, of which 
the magnetic (pyrochlore) sites are the bond midpoints. The 
dipole interaction of spins at r and r' is 



Z? [1 — 3 cos -0 cos -0'] 



„'|3 



(1) 



with D = ^qp? / Anap^ , where p is the moment of one spin, 
and distances are measured in units of a^. Also, cosV" and 
cos -tp' are the angles made by their respective easy axes with 
the vector (r — r'). The effective nearest-neighbor Ising spin 
interaction is Jeff — [J + 5-D)/3 ~ I.IK, with a net ferro- 
magnetic sign that is frustrated since the 3 -fold axes differ 
by the tetrahedral angle 109°, enforcing the well-known two- 
in/two-out spin structure in every tetrahedron. 

I adopt the "dumbbell" approximation for the spin Hamilto- 
nian of the model system; each point dipole is represented by a 
pair of opposite magnetic charges on neighboring tetrahedron 
centers. Then, in any two-in/two-out state, these magnetic 
charges cancel at every diamond lattice site so that (within 
this approximation, and nearly so in dipolar spin ice Ol) all 
such states are exactly degenerate. 

The elementary excitation is a tetrahedron having 3 in and 
1 out-pointing spins, or the reverse, which is the "monopole" 
defect (located at the tetrahedron center). Each monopole con- 
tributes a far field of Coulomb form |B| = (5/47r|R|, where 
Q — ±2/i/arf is its (conserved) magnetic charge, so they are 
indeed emergent mangnetic monopoles f^. Monopole pairs 
interact with the corresponding Coulomb-like potential, and 
even outside the dumbbell approximation each microstate's 
total energy depends almost entirely on the monopole posi- 
tions and very little on the configuration of spins in 2-in/2-out 
tetrahedra apart from the monopoles. where p is the electron 
spin's moment and ad is the separation of tetrahedron centers. 
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A single spin flip in the ground state produces a pair of op- 
positely charged monopoles, which after separating can dif- 
fuse independently over tetrahedron centers or ultimately re- 
combine, much like electrons and holes in a semiconduc- 
tor We let Ao be the cost to create and pull apart a pair of 
monopoles. Within the dumbbell model, this is set by making 
the nearest-neighbor spin interaction to match Jeff, by setting 
Ao = |[J + 4(1 + ^/2jl)D] « 6K li4i]. If we ignore 
Debye screening, the thermally excited density of unbound 
monopoles of both signs is thus 

-Ao/2T 



(2) 



where Vd = v^ap is the volume per diamond vertex, and 
Pauling's approximation for the entropy gives Cmono ~ 4/3 
for the prefactor. 

There are two monopole-dominated temperature regimes. 
Most studies to date (including NMR 1 17]) focused on mod- 
erately low temperatures 1-5 K, at which thermally excited 
monopoles are common but dilute enough to be useful de- 
grees of freedom, but still dense enough that Debye screening 
is crucial and reduces the effective A in (|2]i. The present study 
focuses on the very low temperatures T < T* A/10 ~ 
0.6-ftr, at which rimono < 0.01 and screening is negligible. 
(T* also appears experimentally as the crossover between dif- 
ferent behaviors llisll .) 

Bound pairs of opposite monopoles are possible with fi- 
nite separations. The closest separation is nearest-neighbor 
tetrahedra (separation a^), which is equivalent to a flip of 
just one spin (shared by the tetrahedra) relative to the 2-in/2- 
out background. This has a cost A2 implying a density of 
pairs ~ exp(— A2/T). Although the energy is reduced by the 
Coulomb attraction (A2 < Aq), we lose the factor of 2 in the 
exponent of (|2]i (due to the entropy of deconfinement), so it 
depends on details whether bound or unbound monopoles are 
commoner. As will be shown below, only the unbound ones 
make a novel contribution to the T2 relaxation. (Following 
a quench, "back-to-back" bound pairs of monopoles on adja- 
cent tetrahedra may also persist out of equilibrium, whenever 
the intervening spin has the minority in/out sense. jlSll 

Dy spins are assumed to flip at rate Tq^ if the energy Ai? is 
unchanged or decreases, otherwise t^^ exp(— A£^/T). The 
only spin flip with /S.E — is by one of the 3 spins in 
the majority (in or out) direction on a monopole tetrahedron; 
this transfers the monopole to the other tetrahedron sharing 
the flipped spin. lll9ll Thus the dynamics at low T consists 
of random-walking monopoles which have a combined rate 
to hop of Stq"^, plus the occasional creation or annihilation 
of a monopole pair. Based on magnetic relaxation experi- 
ments flT], it is believed that tq = 10~^ to 10^'^ s and is 
temperature independent. (The tq value can be inferred from 
NMR, as explained further below.) 

NMR model — 

To model the zero-field NMR, I assume a single probe mo- 
ment located at the origin, which for simplicity is imagined to 
be spin 1/2. f20] To zero order, the probe sees a base field ho 
depending on the (frozen) nearby Dy spins, the direction of 



which defines the probe spin's quantization axis. Specifically, 
I have in mind an O^^ nucleus situated at a tetrahedron cen- 
ter: there, ho is due mainly to the four surrounding Dy spins 
(which satisfy the 2-in/2-out constraint) and ho points along a 
(100) axis. The base NMR precession frequency is wo — jho 
where 7 is the gyromagnetic ratio. (In the case of the Ti^'' 
nuclear moment l2li] . wq instead represents its quadrupolar 
frequency.) All nontrivial probe spin dynamics is a response 
to additional, time-varying perturbing fields h(t) due to dis- 
tant flipping spins (presumably Dy: interaction between nu- 
clear moments will be neglected). We call these "flippers" for 
short: they might either be diffusing monopoles, or fixed spin 
impurities, and each flips randomly with a time constant r^. 

Within the material, the frozen local magnetic field has 
strong and complicated spatial modulations within each unit 
cell 022! E3I1 However, only changes of h(t) matter for the 
NMR response: Ti relaxation is sensitive only to high fre- 
quencies — i.e. the flips themselves — while the T2 relax- 
ation is measured using the spin-echo technique which can- 
cels the frozen non-uniformity of the field. Such changes are 
given, at separations R > a^, by the "Coulomb" fields of 
the monopoles, that is ztQ/inR in the radial direction. Ac- 
tually, since a single monopole hop is the same as one spin 
reversal, the field change Ah is just twice the dipole field due 
to that spin, which I describe by a field scale hn/ R^, where 
the distance R is measured in units of Op. We can convert 
the field scale to a frequency using the gyromagnetic ratio, 
cjx = ^J^^djR? ■ The factor \/2 is the root-mean-square of 
the angular factor in ([T]i. 

I will compute both kinds of NMR relaxation functions 
g\{t) and 52(0- To start with, consider the relaxation func- 
tions goi(i) and 502 (^) due to a single flipper at distance |R| 
from the probe spin. To set up the longitudinal or Ti relax- 
ation, say the probe spin is initially aligned with its local axis. 
As the local field direction fluctuates due to the transverse 
components h^ (t) of the perturbing field, the probe spin fol- 
lows adiabatically any slow fluctuations, tilts the field axis 
by a small angle S ss Ah±/hQ; assuming this changes much 
faster than a precession period, the old state is projected onto 
the new axis, so the probability for the probe spin to end up 
with the opposite sense is (1 — cos(5)/2 « (5/2)^. The cor- 
relation of the probe spin with its original sense thus decays 
with time, proportional to the relaxation function 



g,{t) ^ ({cos 6 ff) 



exp 



Ah_L^ t 

Ahl Tx 



(3) 



where the N f is the number of flips in time t, its probability 
given by a Poisson distribution {t/Ty^)^f exp{—t/Tx)/Nfl. 

The "T2" relaxation represents dephasing due to the preces- 
sion of the in-plane angle of tranverse polarization due to the 
fluctuating part of the longitudinal field, uj{t) ~ 7(/i||(i) — 
ho), where 7 is the probe spin's gyromagnetic ratio and 
is field component along the probe spina's local axis. The 
T2 relaxation is measured using the spin-echo technique: the 
polarization is (effectively) flipped by a 7r-pulse at a time t 
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and then evolved till time 2t. The relaxation function is thus 

.92 (2t) = (cos($(2i)) where 



$(2t) = 




is the accumulated phase difference for precession around the 
quantization axis. Thus, the Ti relaxation is dominated by the 
fluctuations at high frequency (compared with ujq), whereas 
the T2 relaxation is dominated by the slower fluctuations at 
frequencies comparable to 1/t. 

Now consider specific ally the dephasing relaxation function 
.902 (R; t) due to a single "flipper"at distance |R|. The change 
in precession frequency at each flip is ±2ajx- If we know that 
exactly one spin flip occurred at a random time within the in- 
terval [0, 2t], it can be worked out that g{2t) = 1 - A{2t), 
where A(2t) = 1 - sin{Lu^t) / (uj^t) « (wxi)^/6. If the flips 
in that time are independent, we use the Poisson distribution 
to obtain 

go2(R;t) = e-^^(*)«e-(*/^-)' (5) 

where Tq^^ — ^uj^t~^ The assumption of independent flips is 
valid for a monopole, ['24]. This is the key difference between 
the spin flip statistics of a monopole and those of a fixed im- 
purity: a conventional Ising impurity must alternate plus and 
minus flips since it has only two states. 

In the case of an alternating flipper, if uj^Tx ^ 1 then 
even one flip by the flipper suffices to randomize the phase, 
so go2(2t) ~ e"'^*/'^>=. On the other hand, if Wx^x ^ 1 and 
t/Tx > 1, then {uj(t')uj(t")) = tj2g-2|t'-t"|/rx and $ has 
a Gaussian distribution with ($(2i)^) = (2t)oj'^Tx, implying 
302 (2t) « e-2*/7'o2 with To-/ = ^x^x. 

We can summarize all cases of the single-flipper results by 

ffo.(R;t)«e-™^»-/l^l^ (6) 

where /3oi = 1, while /3o2 = 1 for a fixed flipper at i > 
Tx, but /?02 = 3 for short times or a diffusing flipper; the 
corresponding time constants are 

Toi^ = {LJoVojofr-'; (7a) 
To2^ = LJD^Tx fixed, Otx; (7b) 
To2^ = {ujof/Tx diffusing or << Tx. (7c) 

and WxTx ^ ojDT'x/IR.p ^ 1 is assumed.' 

Random environments: stretching exponentials due to in- 
homogeneity — 

When different probe moments have different environments 
during the measurement time, the observed signal is an aver- 
age of the relaxation function over these environments and is 
likely to acquire a stretched exponential form. This fact is fa- 
miliar (since the 1960s) in the case that the field fluctuations 
are due to fixed paramagnetic impurities. The same thing hap- 
pens if the fluctuations are due to diffusing monopoles whose 
displacement during the measurement time is small compared 
to their distance |R| from the probe spin. 



I next work through the universal form for the relaxation 
function, averaged over environments, valid for both the Ti 
and T2 relaxations. The key trick, allowing for the simple 
final result, is that (for either kind of relaxation) the combined 
relaxation function due to many flippers at sites {Rj} is sim- 
ply a product g'^°"^^{t) = Y[j 5oj(R; t)- (This follows imme- 
diately from the fact that both kinds of relaxation functions 
are products of independent random variables depending on 
the respective flippers.) A (grand-canonical) ensemble is 
specified by setting the flipper occupation of each site to be 1 
with probability or with probability 1 — rix. 125] Defining 
gi{t) to be the ensemble average of gi°"^^{t), I use conditional 
probability to obtain: 

9^{t) = n ([1 -rix] +«x5o.(R;t)) « e-"''^-^*), (8) 

R 

where, independent of rix 

F,{t) EE ^[i-go.(R;i)]«E[i-^"^'^^'""""$) 

R R 

« ce{t/To^f°^^^ (10) 

where I converted the sum in ^ to an integral, using 
Sr ~^ d^R/A/2 in units of 1/Appl. Phys. ^, thus ce = 
(4V3)/J"(l-e-«')d(l/0 = (327r3/9)i/2 « 10.5. 

We can qualitatively interpret the result ( fTOl i as follows. 
More distant flippers contribute smaller fields which take 
longer to decohere the probe spin; thus in (|9]l, 5oi(R; ~ 
at short |R| or 1 for large |R|. Hence the sum in (|9]l roughly 
counts how many are in the first category; if you write goi ^ 
cxp(— _R^/|R|^), then R* is the effective radius within which 
we count the sites so Fi{t) oc R^. 

Incidentally it can be seen that if Tx is temperature- 
independent, then the only temperature dependence in ^ is 
via rix; thus plots of In gi{t), taken at different temperatures, 
ought to be identical, apart from an overall prefactor which 
gives the temperature dependence of nx(T). 

The final result in all cases is a stretched exponential 

g,{t)^eM-{t/T,f'] (11) 

where (3i — /3oi / 2 and 

T-^' ^ cenxT-/\ (12) 

If the "flippers" (such as monopoles) are thermally excited, 
then the temperature dependence of Ti and T2 follows from 
Ti cx (Tx/nx)^''^v Thus if Tx is temperature independent and 
rix has an activation energy A, the activation energy for Ti is 

A//3.. 

Experimental results and comparison to theory — 
Unpublished work of Kitagawa, Takigawa, et al [T?] found 
stretched exponential forms for both relaxation functions, of 
form (fTTl l The time scales are Ti ^ 1 s and T2 ^ 10^^ s 
at 0.5 K, growing with decreasing temperature. Eventually 
they saturate with Ti ^ 10^ s from 0.2 K downwards, while 
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T2 ~ 10~^ s from 0.4 K downwards. Both exponents (3i de- 
crease with temperature, starting with /3i ~ 1 (i.e. unstretched 
exponential) for T > 1 K, it would appear each exponent satu- 
rates to 1/2 at about the same temperature that the correspond- 
ing Ti saturates, so f3i < P2 at intermediate temperatures. (As 
for the ^'^Ti NMR relaxations, the temperature was not taken 
far below T* so not much can be said about T — > behavior, 
but at the temperatures investigated, both Ti and T2 were an 
order of magnitude longer than they were for ^'^O.) 

The low-T limiting behavior has two fundamental contra- 
dictions with any diffusing monopole theory. Rirst is the tem- 
perature dependence of T2. First, in Eq. (O the only T depen- 
dence comes from the density of flippers (presumed in this 
picture to be monopole density n) or conceivably the flip fre- 
quency T^^. In (fT2l i. this implies Ti must have an activated 
temperature dependence, which contradicts the observed satu- 
ration at low temperatures. (If were also activated, it would 
just add to the activated dependence.) 

The second contradiction is that I found (^2 — 3/2 in the 
monopole diffusion regime, contradicting the experimental 
value ~ 1/2. and indicating that /3o2 ~ 1- But that holds 
only when the flipper is fixed and flipping rapidly compared 
to the measurement time, Tx ^ T2 ~ 10^*s. 

One is forced by the data to assume a density of "flippers" 
that has negligible temperature dependence. This can only be 
some kind of quenched disorder in the material, with a den- 
sity perhaps 10^''-10^"^ per tetrahedron so as to dominate the 
monopole density at temperatures below T*, where the tem- 
perature dependence of Ti and T2 levels off. From here on, 
let us accept that the flippers are fixed, giving /3i = /32 = 1/2, 
and see what the experiment tells us about them. 

From the two relaxation times Ti and T2 we can solve for 
the two unknown parameters Tx and n^. Eliminating from 
Eqs. ©, we find 



My best guess is that extra magnetic (Dy) spins appear on 
the nonmagnetic (Ti) site, as in the "stuffed spin ice" Ii26[ 27 1 



WoT-x ~ V2(Ti/T2)i/2 
-1 - (c^oVT^ 



n 



1/2 



(13a) 
(13b) 



For Dy2Ti207, the two known parameters are ajo/27r « 
20MHz, and wd/wq « 1/20. 

Insertion of these and the experimental T^'s into Eqs. ( fTSl l 
yields Tx ~ 3.5 x 10~^s, which is (barely) consistent with 
the condition Tx <C T2, and w 0.002, still in units of a^^'^. 

What can the fixed flippers be? The out-of-equilibrium, 
back-to-back bound monopole pairs are ruled out: any fluctua- 
tions depend on a momentary energy increase, but the rate ^ 
would be thermally activated, contrary to observation. Dilu- 
tion by nonmagnetic sites produces a very similar situation: 
pairs of unbalanced tetrahedra, each of which is analogous to 
a midgap impurity in a semiconductor The minimum energy 
state adjacent to the nonmagnetic site is a bound pair of half- 
monopoles, and the fluctuation rate is again activated. 

The nearest out-of-plane spins are almost as close, but their 
easy axis is oriented exactly so that the two terms in the dipole 
interaction cancel. 



but much more dilutely. (Note the Ti sites themselves form a 
pyrochlore lattice dual to the Dy pyrochlore lattice.) Indeed, 
it was proposed very recently |28] that in the related mate- 
rial Yb2Ti207, around 5% of the Ti sites are occupied by the 
magnetic Yb ion. 

In order for the "stuffed" spins to fluctuate at such with such 
a short time constant tq, they must be practically decoupled 
from the nearly frozen lattice of regular Dy spins. This de- 
coupling seems plausible when one considers the location of 
the Ti sites, at the centers of hexagons formed by Dy sites, 
and assumes the "stuffed" spin has an easy direction along its 
local three-fold axis. First, the local field at the Ti site due 
to its six nearest-neighbor pyrochlore spins cancels if they are 
all oriented in the same sense around the ring, which should 
happen ^ 12% of the time. Second, the dipole coupling to 
its second-nearest Dy spins happens to have an angular fac- 
tor that exactly cancels. Indeed, it appears from Figure 1 of 
Ref. i22il that the typical local field at a Ti site is ^ 0.2 Tesla 
or about ^ 1 /20 of the maximum local field, which occurs on 
the 0(1) sites containing the probe nuclei. ll29ll 



Conclusion — 

In conclusion, I have rederived the stretched-exponential 
form of the NMR relaxation functions due to independent flip- 
ping spins at random, distant positions, which might be ei- 
ther fixed impurities (weakly coupled to any ordered or frozen 
spin background) or else the spin-flips which induce the hop- 
ping of emergent monopole defects in spin ice. In particu- 
lar, monopole hopping implies a novel power of 3/2 in the 
stretched exponential for the T2 relaxation, in contrast to 1/2 
for a fixed magnetic impurity. 

Experiment at the lowest temperatures shows, for both Ti 
and T2 relaxations, a power tending to 1/2 in the stretched ex- 
ponential and relaxation times tending to a constant, both of 
which are incompatible with the monopole picture. An anal- 
ysis was presented that allows extraction of Tx and from 
Ti and T2 with no bias as to the cause (except it depends cru- 
cially on both kinds of relaxation being due to the same fluctu- 
ations). I suggested that dilute magnetic impurities "stuffed" 
on the Ti sites are responsible. 

It would be interesting to see if NMR relaxation in the 
higher temperature regime around 0.5 K can be explained by 
monopoles. This may be more complex, though: there may be 
no temperature range in which thermally excited monopoles 
are dense enough to dominate over the "stuffed" impurities, 
and are at the same time dilute enough for the small density 
approximations used here. 
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